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Abstract 

Susyno is a Mathematica package dedicated to the computation of the 2-loop renormalisa- 
tion group equations of a supersymmetric model based on any gauge group (the only exception 
being multiple U(l) groups) and for any field content. 

1 Introduction 

Supersymmetry (SUSY) is an elegant and well known extension of the Standard Model (SM) 
which aims at solving or softening some of the SM theoretical and experimental shortcomings. 
Supersymmetric models based on Grand Unified Theories (GUTs) offer the appealing feature of 
being described by a unique gauge group and coupling constant; furthermore, the potential new 
superfields, as well as their dynamics, can provide answers to fundamental problems. For example 
SO(10) based models offer a natural explanation for neutrino masses and mixings via a type-I 
seesaw. This occurs since the spinor representation of the group contains a SM singlet superfield 
(with the appropriate quantum numbers). Even larger and/or more complex groups are often 
considered, for example Eq in string inspired models. 

The analysis of the theoretical and phenomenological implications of SUSY GUT models re- 
quires a careful study of the evolution of the fundamental parameters from the high-energy scale 
down to the electroweak scale, at which observables are computed and constraints applied. As 
such, knowledge of the renormalisation group equations (RGEs) is necessary. Although the RGEs 
of several models (e.g., MSSM, NMSSM) are already known [HE], for other SUSY extensions of 
the SM complicated general equations must be used [U[2]. 

Here we describe Susyno, a Mathematica-based package that addresses this issue. The program 
takes as input thegauge group, the representations (i.e., the chiral superfield content), the number 
of flavours/copies^] of each representation, and any abelian discrete symmetries (e.g., R-parity). 
Susyno then constructs the most general superpotential and soft SUSY breaking Lagrangian con- 
sistent with the field content and symmetries imposed. Once these elements have been derived, 
Susyno then computes the 2-loop /3-functions of all the parameters of the model, which is its main 
output. 

throughout this document, the word flavours is meant to convey the sense of copies or repetitions of a given 
representation /field . 



There is another Mathematica package, SARAH [4j, which works in a similar way for models 
based on SU(n) gauge factor groups. However, Susyno is prepared to accept as working input any 
gauge group that does not contain more that one U(l) factoiQ 

This document is organised as follows. Section [2] explains how to install the program and run 
a first simple example (MSSM based). Sections [3] and [4] explain how to prepare the input and 
how to read and interpret the output, also using as examples the MSSM case. Finally, Section [5] 
summarises the tests conducted to validate the code. 



2 Installation and quick start 

Susyno works on Windows, Linux and Mac OS provided that Mathematica 7 (or a latter version) 
is installed. The program is obtainable from 



http : / /web . ist . utl . pt/renato . f onseca/ susyno . html 



The files LieGroups.m, SusyRGEs.m and Susyno. m are the core of the program. These and 
other auxiliary files can be found inside the folder Susyno, which must be extracted from the 
downloaded Susyno- 1.0. tar file to a location that is visible to Mathematica. Typing $Path in 
Mathematica will show a complete list of acceptable locations. A good choice is to place the whole 
folder (not just its contents) in 

(Mathematica base directory) / 'AddOns/ 'Applications 

(note that in a Windows system the slashes "/" must be replaced by backslashes "\"). The package 
can be loaded by typing 



« Susyno" 



in Mathematica's front end. A text message is returned, informing that a built-in help system 
provides a detailed description of the program and its functions (see Appendix A). A tutorial is 
also included. 

The Susyno lines below allow a simple and easy first run: the example consists in a possible 
way of writing the MSSM input. 



groups = {SU2 , SU3}; 



normalization = Sqrt [3/5] ; 

Q = {1/6 normalization, {1}, {1, 0}}; 

u = {-2/3 normalization, {0}, {0, 1}}; 

d = {1/3 normalization, {0>, {0, 1}}; 

L = {-1/2 normalization, {1}, {0, 0}>; 

e = {normalization, {0}, {0, 0}}; 

Hu = {1/2 normalization, {1}, {0, 0}}; 

Hd = {-1/2 normalization, {1>, {0, 0}}; 

2 Multiple U(l) groups would lead to "U(l) mixing" and at present there are still no complete published results 
for this case [5|. 
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content = {Q, u, d, L, e, Hu, Hd}; 



nFlavours = {3, 3, 3, 3, 3, 1, 1}; 
discreteSym = {-1, -1, -1, -1, -1, 1, 1}; 

BetaFunctions2L [groups , content, nFlavours, discreteSym] 

Evaluation of this simple code generates the 2-loop f3- functions of the model (MSSM in this case). 
Notice that no external input or output files are used - everything happens on Mathematica's front 
end. 

3 Defining a model 

A SUSY model contains two building blocks: a superpotential and a soft SUSY breaking La- 
grangian. For a simple gauge group, a general superpotential can be written as 

W = ±Y iik $ i $ j $ k + ^//'"l>,<l> ; + L% , (1) 

with <]? denoting a chiral superfield, and where Y, \i and L are dimensionless, mass, and mass 2 
parameters, respectively. A generic soft SUSY breaking Lagrangian reads 




where A is a gaugino field and (j) denotes scalar fields. M (m 2 ) correspond to gaugino (scalar) soft 
breaking masses, while h, b and s denote trilinear, bilinear and linear soft breaking parameters. 
Susyno works as follows: it requires as input the gauge group, the representations/fields, the num- 
ber of flavours of each representation/field and the discrete abelian symmetries (if any is present) 
of the model. With this information the program then builds internally the superpotential and 
the soft SUSY breaking Lagrangian using an algorithm to automatically name the parameters of 
the model (see Section 4). 

Assuming that this information is summarised in groups, content, nFlavours and discreteSym, 
the user must then call the function BetaFunctions2L (or BetaFunctionslL for 1-loop results) as 
follows: 

BetaFunctions2L [groups , content, nFlavours, discreteSym] 

We shall focus first on each of these arguments. The syntax of the functions will be discussed 
at the end of this section. 

3.1 Gauge group 

The program needs a complete list of all the simple Lie groups of the modeH. For the MSSM (for 
completeness, the superpotential and SUSY soft breaking Lagrangian of the MSSM are presented 
in |Appendix B[ ) this would correspond to 



3 We emphasise here that we are actually dealing with algebras, not groups. Nevertheless, we will adopt the 
common practice in high-energy physics of using group for both these concepts. 
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groups = {SU2, SU3}; 



Any simple group can be given as a factor: the simple gauge factor groups, as well as their 
corresponding Susyno input are collected in Table [TJ 

Simple gauge factor group Susyno input 



SU(n) SU2, SU3, SU4, SU5, ... 

SO(rc) S03, SQ5, S06, S07, ... 
Sp(2n) SP2, SP4, SP6, ... 

G 2 G2 
F 4 F4 
Eg, E7, Eg E6, E7, E8 

Table 1: Simple gauge factor groups 

Notice that despite its generality, Susyno only assigns to variables a finite number of groups 
(for instance variable S0100 is not set). Should this ever become a problem, there is an easy 



method of circumventing it, discussed in Appendix C 



Since it is not a simple group, U(l) is a special case and should not be included in the list of 
simple factor groups. The information regarding any U(l) factors present in the gauge group of 
the model is given through the number of hypercharges of the fields (see below) . 



3.2 Representations/fields: the content of the model 

As mentioned before, Susyno is designed to accept an arbitrary field content. An input must be 
provided in the form of a list 



content={f ieldl ,f ield2, . . .}; 



containing all the gauge group irreducible representations present in the model. It is important to 
stress that the ordering of the fields in content is arbitrary, and decided by the user. However, 
once content has been read by Susyno, the user must consistently adhere to the chosen ordering 
when inputting lists composed of field attributes (for example the number of flavours, as will be 
discussed in the following subsection). Each field should be specified as 

f ield={hChrgl ,hChrg2, . . . ,hChrgM,rep_simplegroupl,rep_simplegroup2, . . 

The first entries correspond to the hypercharges of field (if any)@, which are just real numbers. 
After the hypercharges one must declare the representations of field under each of the simple 
gauge factor groups mentioned above. These representations must be specified by their Dynkin 
coefficients (see |Appendix C] for details). In Table [2] we list some of the representations of SU(2), 
SU(3), SU(5) and SO(10) (Dynkin coefficients and corresponding dimensions). 

To understand how the MSSM was specified in the example of Section [21 we just need the 
following information from Table [2 

• the Dynkin coefficients of the trivial and fundamental representations of SU(2): {0} and {1}; 

4 Notice that, while it is possible to have more than one hypercharge (i.e. more than one U(l) factor group), the 
formalism of Susyno is based on [UH], which do not include the effects of a possible U(l) mixing. 
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Group Representation 

SU(2) 

SU(3) 
SU(5) 



SO(10) 



Dynkin coefficients 


Dimension (Name) 


10} 


1 (Trivial/Singlet) 


{1} 


2 (Fundamental/Doublet) 


L J 


3 (Adjoint/Triplet) 


{0,0} 


1 (Trivial/Singlet) 


{1,0} 


3 (Fundamental) 


{0,1} 


3 (Anti-fundamental) 


{1,1} 


8 (Adjoint) 


{0,0,0,0} 


1 (Trivial/Singlet) 


{1,0,0,0} 


5 (Fundamental) 


{0,0,0,1} 


5 (Anti-fundamental) 


{0,1,0,0} 


10 


{2,0,0,0} 

L ' ' ' J 


15 


{0,0,0,2} 


15 


{1,0,0,1} 


24 (Adioint) 


{0,0,0,0,0} 


1 (Trivial/Singlet) 


{1,0,0,0,0} 


10 (Fundamental) 


{0,0,0,0,1} 


16 (Spinor) 


{0,0,0,1,0} 


16 (Spinor's conjugate) 


{0,1,0,0,0} 


45 (Adjoint) 


{2,0,0,0,0} 


54 


{0,0,1,0,0} 


120 


{0,0,0,0,2} 


126 


{0,0,0,2,0} 


126 



Table 2: List of some of the representations of SU(2), SU(3), SU(5) and SO(10) 



• the Dynkin coefficients of the trivial, fundamental (3) and anti-fundamental (3) representa- 
tions of SU(3): {0,0}, {1,0}, {0,1}. 

For the MSSM each field must then be cast in the format 



f ield={Ul_charge , SU(2) _rep , SU (3) _rep} ; 

Further normalising the hypercharges with the usual y^f factor (from an embedding of the MSSM 
in an SU(5) based model), we can then write the following: 



normalization Sqrt[3/5]; 

Q = {1/6 normalization, {1}, {1, 0}}; 

u = {-2/3 normalization, {0}, {0, 1}>; 

d = {1/3 normalization, {0}, {0, 1}}; 

L = {-1/2 normalization, {1}, {0, 0}>; 

e = {normalization, {0}, {0, 0}}; 

Hu = {1/2 normalization, {1}, {0, 0}>; 

Hd = {-1/2 normalization, {1>, {0, 0}}; 

content = {Q, u, d, L, e, Hu, Hd}; 
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We must emphasise here that although the user is free to choose the ordering of the simple 
factor groups, SU(2) and SU(3), please bear in mind that once this is set (e.g., groups = {SU2, 
SU3}) one must adhere to the (user-established) convention, and define the representations of the 
fields accordingly: 



f ield={Ul_charge , SU(2) _rep , SU (3) _rep} ; 



Also notice that Susyno accepts any choice for the normalisation of the hypercharges. 



3.3 Number of flavours and abelian discrete symmetries 

Susyno needs two more input lists: one containing the number of flavours of each field and another 
defining its abelian discrete symmetries. As previously mentioned, the attributes of the fields listed 
in both nFlavours and discreteSym must follow the field ordering established by content. In 
our MSSM example we used the ordering {Q, u, d, L, e, Hu, Hd}, and therefore 

nFlavours = {3, 3, 3, 3, 3, 1, 1}; 
discreteSym = {-1, -1, -1, -1, -1, 1, 1}; 

In this particular case, it is clear that the discrete symmetry imposed corresponds to R-parity. For 
the most general (R-parity violating) MSSM we have 



discreteSym = {1, 1, 1, 1, 1, 1, 1>; 



Let us consider another example: for instance, if we were to modify the MSSM to include m 
copies of H u and H^, we would write 



nFlavours = {3, 3, 3, 3, 3, m, m}; 



Both these lists can be empty: if the flavour list is {}, the program will assign by default the 
symbolic values {nf [1] ,nf [2] ,...}; if the discrete symmetry list is {}, Susyno will consider it to 
be {1 , 1 , . . . }, meaning it will assume that all fields transform trivially. 

3.4 Calling the functions BetaFunctionslL, BetaFunctions2L 

The four lists discussed so far are the arguments of BetaFunctionslL or BetaFunctions2L, 

BetaFunctionslL [groups , content, nFlavours, discreteSym, <print>] 
BetaFunctions2L [groups , content, nFlavours, discreteSym, <print>] 



These functions compute the one- and two-loop RGEs of the model (see Appendix D). The op- 
tional argument <print> can be set to True (default value) or False. If <print>=True the results 
are printed on the screen but the functions themselves do not return a value. On the contrary, if 
<print>=False there is no printing and the functions now return the results as a list. Assuming 
that the model contains parameters named XI, X2, the format of this output list is 

{{XI, /3^ ) >,{X2, P$y,...} (BetaFunctionslL) 
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{{XI, 4i ) >,{X2 ) Z?^, ^2 ) >»---> (BetaFunctions2L) 

that is, a set of sublists each containing a parameter Xi, and the one-loop (one- and two- loop) 
beta functions for its RG evolution: (3^ (/3$ , 

If for some reason the results are to be further processed with Mathematica, then the best 
choice is clearly to set <print>=False. 

4 The output of Susyno 

Once all the input is provided, Susyno automatically builds the Lagrangian of a model, neither 
asking for "names" of fields nor "names" of parameters. On the one hand, inputting a model 
becomes very easy, since it is not even necessary to know the exact number of its parameters. On 
the other hand, this renders the output hard to read (and hence not particularly user- friendly) , 
since the names of the parameters are chosen by the program. 

Still, provided we know the conventions used by the program, the notation used by Susyno can 
be manually changed by the user (after the code is executed). We notice that in the present version 
of Susyno there are no custom built-in functions to export the results. Users who wish to do so 
must do it manually or with the help of Mathematica's built-in functions CForm and FortranForm. 

4.1 Naming of parameters 

Susyno assigns names to the parameters of a model in such a way that the user can identify which 
representations/fields they are multiplying: 

y[fieldl, field2, field3, Invlndex, <flavl>, <flav2>, <flav3>] 
mu[fieldl, field2, Invlndex, <flavl>, <flav2>] 
l[fieldl, Invlndex, <flavl>] 

h[fieldl, field2, field3, Invlndex, <flavl>, <flav2>, <flav3>] 
b[fieldl, field2, Invlndex, <flavl>, <flav2>] 
s[fieldl, Invlndex, <flavl>] 

m2[fieldl, field2, Invlndex, <flavl>, <flav2>] 

A few comments concerning the above (output) tensors are in order: 

• y, mu, 1, h, b, s and m2 can be easily identified with the different types of couplings and 
dimensionful parameters of the superpotential and the soft SUSY breaking Lagrangian; 

• f ieldl, f ield2, f ield3 are the indices of the fields entering a given coupling, or soft breaking 
term. For example, with content={Q, u, d, L, e, Hu, Hd} we have Q=l, u=2, d=3, L=4, 
e=5, Hu=6, Hd=7. The up-quark Yukawa couplings would then be y [1,2,6, ...] ; 

• There is the possibility that the product of 3 representations, Ri <g> R2 <g> R3, contains more 
than one invariant. Therefore we need Invlndex=l,2, ... to distinguish them. This is rare 
though, so in most cases (e.g., the MSSM) Invlndex=l for all parameters. Notice that in 
linear and bilinear terms, R\ and R\ <g> i?2, this problem does not arise since there is at most 
one invariant. However, for consistency Susyno still uses Invlndex in this case, although 
setting its value to 1; 
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• <flavl>, <flav2> , <flav3> are the flavour indices of fieldl, field2, field3. If any of 
these fields has only one flavour, the corresponding index is omitted. Consider again the 
example of the up-quark Yukawa couplings: we would have y[l,2,6,l,i,j,k] where i 
flavour of Q, j = flavour of u, k = flavour of Hu. Yet Hu only has one flavour so the correct 
parameter name is y [1 , 2 , 6 , 1 , i , j] . 

Additionally, there are also the coupling constants and the gaugino masses: 

g[l], g[2], ... 
M[l] , M[2] , . . . 

A list containing the matching of MSSM's standard notation parameters and Susyno assigne- 
ments can be found in |Appendix B| 

4.2 The full Lagrangian as used by the program 

Knowing the full Lagrangian can be important. As an example, consider a possible normalisation 
issue: how can we be sure that Susyno is correctly taking fJ,H u ■ as mu [6 , 7 , 1] H u ■ H^, and 
not -mu [6 , 7 , 1] H u ■ H^, 2mu [6 , 7 , 1] H u ■ or any other multiple? The only way to solve such 
problems related to the normalisation of the parameters is to check the Lagrangian. This can be 
done using the ShowLagrangian function: 

ShowLagrangian [groups , content, nFlavours, discreteSym, <print>] 

The first four arguments are the same as for the BetaFunctionslL and BetaFunctions2L 
functions. As for the last one, setting <print>=True (default value) will print the different parts 
of W and ^f so ft- 

• ^•'• /, <lvlvl>, (Y part), 5// ; "I>,<l> ; (/i part), L% (L part); 

• \h^ k ^i(\)j(\) k (H part), \W fafa (B part), s^i (S part) and (m 2 ) (M2 part). 

If <print>=False no printing is done; instead the function returns the list of tensors {Y, fx, L, h, b, 
s,m 2 |. Note that Susyno never expands the flavour indices of these tensors so the first, second 
and third indices are assumed to carry symbolic flavours i (or none), j (or none), k (or none) 0. 

ShowLagrangian is completely independent of BetaFunctionslL and/or BetaFunctions2L, 
and so it can be called at any time to generate the superpotential and soft SUSY breaking La- 
grangian of a given model. 

Finally, notice that in order to build Jz? it is necessary to assume a particular basis for the 
representations/fields of the model. Quantities such as the (3- functions clearly do not depend on 
this choice. However, it is important to know that Susyno can use unconventional basis for the 
representations. The following is a good example. According to the program, two fields, A={a[l], 
a[2], a[3]} and B={b[l], b[2], b[3]}, transforming as 3 and 3 of SU(3) will form the invariant 
a[l]b[3] - a[2] b[2] + a[3]b[l], not a[l]b[l] + a[2] b[2] + a[3]b[3]. Both are valid expressions as 

5 These flavour indices do not show up if the particular entry of the tensors we are considering is unflavoured. For 
example the entry of the Y tensor related to the fields Q, d, Hd will depend on i (flavour of the first field, Q) and 
j (flavour of the second field, d) but not on k since the third field (Hd) is unflavoured. Notice also that exchanging 
for example Hd with Q corresponds to having a different entry of the Y tensor: this is because the flavour of Q 
(now the third field) is assumed to be symbolically k instead of i. 
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long as they are used in a consistent way - it is just a matter of the choice of basis. Thus, when 
analysing the output of ShowLagrangian we must distinguish between irrelevant basis choices and 
relevant normalisation conventions for the parameters. The /3-functions are only sensitive to the 
latter ones. 

5 Tests/validation of Susyno 

The output of Susyno was confronted with the analysis of some models available in the literature. 
In particular, the RGEs generated by Susyno were compared with the results of [T] (MSSM), [5] 
(R-parity violating MSSM), [3] (general NMSSM) as well as [7] (SU(5)-based models). 

In general the program's RGEs are consistent with the results collected in the above publi- 
cations, but in some cases differences were found. These have been collected in the program's 
webpage. 
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Appendix A List of available functions 

As in most programs, the code Susyno is spread over many internal functions. Due to their nature, 
some of these functions may be useful on their own, and they were thus built in a user-friendly 
way, and are documented. 

Below is the list of functions that can be called directly by the user in Mathematica's front-end, 
followed by a brief description. The built-in help system describes in detail how to use them. 

• Adjoint - Computes the Dynkin coefficients of the adjoint representation of a group; 

• BetaFunctionslL - Computes the 1-loop /3-functions of a SUSY model; 

• BetaFunctions2L - Computes the 2-loop /3-functions of a SUSY model; 

• CanonicalForm - Simplifies an expression written in Einstein's notation; 

• Casimir - Computes the quadratic Casimir of a representation; 

• CM - Computes the Cartan matrix of any simple group; 

• DimR - Computes the dimension of a representation; 

• Invariants - Computes (in a given basis) the invariant combination(s) of a product of one, 
two or three representations; 

• ListContract - Efficiently calculates traces of multi-index sparse tensors; 

• PositiveRoots - Computes the positive roots of a group; 

• ReduceRepProduct - Decomposes a direct product representation in its irreducible parts [8]; 

• RepMatrices - Computes (in a given basis) the explicit representation matrices; 

• ShowLagrangian - Generates the Lagrangian [W and ^£ S ojt) of a SUSY model; 

• Weights - Computes the weights of a representation, including degeneracy. 

Appendix B Parameters of the MSSM 

For completeness, we include here the superpotential and soft SUSY breaking Lagrangian of the 
MSSM, following the conventions of [lj. 



In the above, and as usual, hats denote superfields and tildes the scalar component of each chiral 
superfield. The contraction of two SU(2) doublets (denoted by •) follows the convention H u ■ Hd = 
H^Hj — H®Hj. In addition to the parameters entering W and Jt? S oft, one still has the gauge 
couplings g a . Table Q] summarises the matching of the physical notation of Eqs. ([3|4p with the one 
used by Susyno. 




W 



(3) 



(4) 
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Parameter Susyno's notation 

91, 92, 93 g[l], g[2], g[3] 

Mi, M 2 , M 3 M[l] , M[2] , M[3] 
y[l,2,6,l,j,i] 

(YdJy y[l,3,7,l,j,i] 

(F e ) . y[4,5,7,l,j,i] 

/i mu[6,7,l] 

(hu)^ h[l,2,6,l,j,i] 
h[l,3,7,l,j,i] 

(he)^ h[4,5,7,l,j,i] 

B b[6,7,l] 

(m 2 Q y m2[l,l,l,j,i] 

(m 2 u ).. m2[2,2,l,i,j] 

(ml).. m2[3,3,l,i,j] 

(ml).. m2[4,4,l,j,i] 

(mfj^ m2[5,5,l,i, j] 

m 2 H m2 [6,6,1] 

m 2 rr m2 [7,7,1] 

Table 1: Parameters of the MSSM assuming the field ordering {Q,u,d,L,e,Hu,Hd} and the simple 
factor group ordering {SU2,SU3}. 



Appendix C More on groups and representations 

A simple complex Lie algebra is specified by a Cartan matrix [9] which is a square matrix with 
diagonal entries all equal to 2 and off-diagonal entries equal to 0, -1, -2 or -3. For example, the 
Cartan matrix of SO (10) is 





2 


-1 











\ 




-1 


2 


-1 















-1 


2 


-1 


-1 












-1 


2 







V 








-1 





2 


/ 



The Cartan matrix is the ingredient required by Susyno whenever a group is being asked as input. 
Although the user can manually provide the Cartan matrices, the program already has a built-in 
function, CM, that can compute them for any simple group. The syntax is 

CM["group name",f amilylndex] 

where "group name" can be "SO", "SU", "SP", "G", "F" or "E" and familylndex should be an 
integer. For example 

CM ["SO", 10] 

will return the above 5x5 Cartan matrix. This function will work for virtually any simple 
group. However, this is still a heavy syntax, so the program assigns to some variables (SU2, SU3, 
SU5, S010, etc.) the corresponding Cartan matrices (e.g., S010 is the same as CM["S0", 10]). For 
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the case of SO(ra), SU(n) with n > 32 and Sp(2m) with m > 16 no variables were set, and the CM 
function should be used. 



As mentioned in the text, Susyno requires that each field be specified as a list of hypercharges 
and representations under the simple gauge factor groups. We will focus now on this last element: 
representations of simple groups. 

A representation of a (simple) group is often labeled by its dimension (e.g., 3 of SU(3), 54 of 
SO(10), ...) but this can be ambiguous. To avoid the problem, we should use Dynkin coefficients 
which are a list of n non-negative integers. This number n is known as the group rank and is equal 
to the size of the (square) Cartan matrix of the group. For instance, since the Cartan matrix of 
SO(10) is 5 x 5, all representations of the group are of the form 



{il,i2,i3,i4,i5} 



(il,.„ 5 i5 are non-negative integers). 

There are functions in Susyno that compute properties of the representations (e.g., DimR cal- 
culates the dimension of a representation, ReduceRepProduct reduces products of representations) 
and they are documented in the built-in help files (see also Appendix A ) . These should be enough 
to allow recognising a representation by its Dynkin coefficients; however, should the user wish 
to consult lists of representations, these are available in the literature [10J. Some of the most 
frequently used representations were shown in table [2j 



Appendix D Renormalisation group equations of a general model 



Susyno uses the generic two loop RGEs of [TJ [2] . For completeness these results are reproduced 

here (for a model based on a simple gauge group G). 

In general a parameter X evolves according to the equation 

dt 16ir 2 ^ x (16^ 2 ) 2 

where t = logQ (Q being the renormalisation scale parameter). The /3-functions for the sev- 
eral parameters (gauge coupling constants, superpotential parameters and SUSY soft breaking 
parameters) are as follows. 

Coupling constants 

UP = g 3 [S (R) - 3C (G)] (7) 
4 2) = 9 5 {-G[C(G)} 2 + 2C(G)S(R) + 4S(R)C(R)}-g 3 ^Y^ k Y l , k (8) 
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Superpotential parameters 



where 



ijfe 



r^ 7 p + (fe o ») + (fe o j) 



v( 2 )i 



pf 



= /i*7^' + (J <) 
= 



7 j 1}< = -^ M y^-2^'c(») 

7f ]i = -\Y imn Y n ™Y pqr Y mr i + S 2 y iM P' [2G (p) -C(i)} 
+25jg 4 \c (i) S (R) + 2C (i) 2 - 3G (G) C (i) 



Soft SUSY breaking Lagrangian parameters 



Pm 

o(2) 

Pm 



g 2 [2S (i?) - 6C (G)] M 



-24G (G) 2 + 8G (G) 5 (i?) + 165 (R) C (R) 



(9) 
(10) 

(11) 
(12) 

(13) 
(14) 

(15) 
(16) 



M + 2g- 



C(k) 
d(Gj 



(17) 

(h ijk - MY ijk ^j Y ijk 
(18) 



ft 



(i) 



ft 



(2) 



ijk 



ijk 



h^ l Y lmn Y mnk + y^ mn /i mnfe - 2 - 2Mr^'j g 2 c(k) + (k^i) + (k^ j) 

(19) 



— hV^V, V n PQV vmrk_-yijl^ yripq-y yrark _ vijlxr Lnpqy -ymrk 
r,' 1 1 lmn 1 1 pqr 1 1 1 lmn 1 1 pqr 11 ' 1 1 lmn 11 ' 1 pqr 1 



+ 
+ 



(h ijl Y lpq YPi k + 2Y^ l Y lpq h pqk - 2MY i i l Y lpq Y™ k } g 2 [2C(p) - C(fc)] 

(2/1^ - 8Mr^) 5 4 [G(fc)S(i?) + 2C(£;) 2 - 3C(G)C(k)] + (jfe <-> i) + (jfe <-> 7) 



(20) 
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(1) 



■IJ 



h il Y lmn Y™i + \ Y ^Y lmn b mn + fi il Y lmn h m ^ - 2 (&« - 2M^) g 2 C(i) + (i« j) 



(21) 



(2) 



— -h il Y, V pqn Y Y mr i — -Y^ l Y h mr V Y pqn — -Y^ l Y ii mr Y h pqn 
2" 1 lmn 1 1 pgr 1 ^ lmn u 1 pqr 1 ^ 1 lmnf J ' 1 pqr' L 

-fi U Y lmn h npq Y pqr Y m ^ - v il Y lmn Y npq Y pqr h mr i + 2Y^ l Y lpq (&w - fi pq M) g 2 C(p) 

+ (V% g y^ + 2//%,/^ - 2//% ? y^M) g 2 [2C( P ) - c(i)] 

+ (2b ij - 8n ij M) g 4 [C(i)S(R) + 2C(i) 2 - 3C(G)C(i)] + (i <-> j) (22) 



-y ip % grS r + /^y^z/ + + 2Y imn { m 2 ) l M + m pg (23) 

2<? 2 C (q) Y^Y pqr s r - ^Y irn Y npq Y™ m Y mrl s l - 4g 2 C (q) (Y ipq M - h™) Y pql L l 

- {Y mn Y npq h™ r Y rml + h imn Y npq Y^Y rml ) L l - Ag 2 C (q) tf l Y lpq (fi pq M - 

-fi u (Y lmn h npq Y pqrf i rm + Y lmn Y npq Y pqr b rm ) + Ag 2 C (q) [2Y^^ pq \M\ 2 - Y ipq b pq M 

-h^v pq M* + h^b pq + Y^(m 2 Y pt i rq + Y^{m 2 )lii pq ] - Y imn Y npq hP% m 



uimny Y pql hi _Y imn h h^ 1 1 



h imn h nnn Y pql i 



' n\«j > Olm * n>npqH' ' f^lm ^ ,lj npq^ ' ' Him Y {jll ^ m IJ,i r Y Pq Y pqn 

-Y* mn Y npq Y pql (m 2 y if i rm - Y* mn (m 2 ) l n Y lpq Y pqr v rm - 2Y* mn Y npq (m 2 ) q Y rpl fi lm 

(24) 



l -Y ipq Y^ n {m 2 ) j n + l -Y^Y pqn {m 2 ) r : + 2Y ipq Y^{m 2 )l + h ipq V™ 



-8% \M\ 2 g 2 C(i) + 2g 2 tfTi (t A m 2 ) 



(25) 



— -(m 2 \ l Vi Y mr iY Y pqn — -(m 2 yY lmn Y Y pqr Y —Y, Y^ nm (m 2 \ l Y Y rpq 

<2^ )i 1 lmn 1 1 pqr 1 2 1 mri 1 1 pqn 1 ilm 1 V'" 1 )r n P1 

_V-, Yi nm (m 2 \ r Y Y lpq — V, Vi nr (rn 2 \ l V VPQ m OYi vi ln V Y mpr (m 2 \ q 

1 ilm 1 \i' b )n r P1 1 tlm 1 V'" 1 )n VQ r ^ 1 ilm 1 1 npq 1 y" 1 ) r 

—Y, Yi ln h h mpq — h-, W ln Y Y mpq — h-, Yi ln Y h mpq — Y, W ln h Y mpq 

1 ilm 1 "'npq 11 ' 'Him 1 '' 1 npq 1 'Him 1 1 npq"' 1 ilm' 1 ' ""ripq 1 

+ {m 2 )\Y lpq Y^ + Y ipq Y l P q (m 2 )i + AY ipq Y^ l {m 2 ) q + 2h ipq h™ q - 2h ipq Y^M 

-2Y ipq h^M* + 4Y ipq Yi™ \M\ 2 g 2 [C(p) + C(q) - C(i)] - 2g 2 tf\t A m 2 ) l r Y lpq Y^ 

+8/tf j Tr [t A C{r)m 2 ] + Sjg 4 \M\ 2 [24C(i)S(R) + 48C(i) 2 - 72C{G)C{i)] 
+85lg 4 C(i) {Tr [S(r)m 2 ] - C{G) |M| 2 } (26) 



A few comments and clarifications are still in order concerning the variables appearing in the 
different /3-functions: 

. Y ijk = (Y^y, h ijk = (h^y, m 6 = (v ij y and bn = {b^y- 

• d (G) — Dimension of the adjoint representation of group G; 

• C (i) — Quadratic Casimir invariant of the representation of the chiral superfield with index 
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• C (G) — Quadratic Casimir invariant of the adjoint representation of group G; 

• S (R) = Dynkin index summed over all chiral multiplets. However S (R) C (R) should be 
interpreted as the sum of Dynkin indices weighted by the quadratic Casimir invariant; 

• t A — Representation matrices under the gauge group G. Terms with t A are only relevant 
for U(l) groups; 

• In B I and B^ 2 \, the traces should be understood as traces over all chiral superfields. 

If the gauge group is not simple but rather a direct product of simple groups and U(l) factors, the 
equations shown here must be adapted. The modifications are fairly straightforward, and can be 
found in [lj. 
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